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Abstract — We investigate the problem of design- 
ing a minimum cost flow network interconnecting n 
sources and a single sink, each with known locations 
and flows. The network may contain other unpre- 
scribed nodes, known as Steiner points. For concave 
increasing cost functions, a minimum cost network of 
this sort has a tree topology, and hence can be called a 
Minimum Gilbert Arborescence (MGA). We charac- 
terise the local topological structure of Steiner points 
in MGAs for linear cost functions. This problem has 
applications to the design of drains, gas pipelines and 
underground mine access. 

Keywords: optimisation, networks, network flows, 
Steiner trees 

1 Introduction 

One of the most important advances in physical network 
design optimisation since the 1960's has been the devel- 
opment of theory for solving the Steiner Minimum Tree 
(SMT) problem. This problem asks for a shortest net- 
work spanning a given set of points, called terminals^ in 
a given metric space. It differs from the minimum span- 
ning tree problem in that additional points, called Steiner 
•points, can be included to create a spanning network that 
is shorter than would otherwise be possible. This has nu- 
merous applications, including the design of telecommu- 
nications or transport networks for the problem in the 
Euclidean plane (the I2 metric), and the design of mi- 
crochips for the problem in the rectilinear plane (the li 
metric) [5j. 

Gilbert W proposed a generalisation of the SMT prob- 
lem whereby symmetric non-negative flows are assigned 
between each pair of terminals. The aim is to find a least 
cost network interconnecting the terminals, where each 
edge has an associated total such that the flow condi- 
tions between terminals are satisfied, and Steiner points 
satisfy Kirchhoff's rule (ie, the net incoming and outgo- 
ing flows at each Steiner point are equal). The cost of 
an edge is its length multiplied by a non-negative weight. 
The weight is determined by a given function of the total 
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flow being routed through that edge, where the function 
satisfies a number of conditions. The Gilbert network 
problem (GNP) asks for a minimum-cost network span- 
ning a given set of terminals with given flow demands and 
a given weight function. 

An important variation on this problem that we will show 
to be a special case of the GNP occurs when the terminals 
consist of n sources and a unique sink, and all flows not 
between a source and the sink are zero. This problem has 
applications to drainage networks [7], gas pipelines [1], 
and underground mining networks [2]. 

If the weight function is concave and increasing, the re- 
sulting minimum network has a tree topology, and pro- 
vides a directed path from each source to the sink. Such 
a network can be called an arborescence, and we refer 
to this special case of the GNP as the Gilbert arbores- 
cence problem (GAP). Traditionally, the term 'arbores- 
cence' has been used to describe a rooted tree providing 
directed paths from the unique root (source) to a given 
set of sinks. Here we are interested in the case where the 
flow directions are reversed, i.e. flow is from n sources 
to a unique sink. It is clear, however, that the resulting 
weights for the two problems are equivalent, hence we 
will continue to use the term 'arborescence' for the latter 
case. Moreover, if we take the sum of these two cases, and 
rescale the flows (dividing flows in each direction by 2), 
then again the weights for the total flow on each edge are 
the same as in the previous two cases. This justifies our 
claim that the GAP can be treated as a special case of 
the GNP. It will be convenient, however, for the remain- 
der of this paper to think of arborescences as networks 
with a unique sink. 

A minimum Gilbert arborescence (MGA) is a (global) 
minimum-cost arborescence for a given set of terminals 
and flows, and a given cost function. All flows in the net- 
work are directed towards the unique sink. MGAs have 
been used to model drainage networks, such as plumb- 
ing networks in buildings [7l , as well as gas pipeline net- 
works pi. In both cases heuristics were provided for ob- 
taining low-cost solutions. A special complication for 
these networks is that the cost of an edge depends on 
the diameter of the pipe, and the diameter of the pipe 
depends on its length. This complication will not be con- 
sidered in this work. 

More recently, MGAs have been used to model under- 
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ground mining networks [5] . Given a set of underground 
locations, called draw points, and their estimated ore 
tonnages, the development and haulage costs associated 
with an underground mine can be minimised by finding 
an MGA interconnecting the underground points with 
a fixed breakout point at the surface. Although the Eu- 
clidean metric is useful in some situations, more generally 
additional constraints must be imposed for truck naviga- 
bility such as a gradient constraint, which can be handled 
by altering the metric |2j- 

In Section 2 we specify the nature of the weight func- 
tion that we consider in this paper, and formally define 
minimum Gilbert networks and Gilbert arborescences in 
Minkowski spaces (which generalise Euclidean spaces). In 
Section 3 we give a general topological characterisation 
of Steiner points in such networks, for smooth Minkowski 
spaces. We then apply this characterisation, in Section 4, 
to the smooth Minkowski plane with a linear weight func- 
tion to show that in this case all Steiner points have de- 
gree 3. 

2 Preliminaries 

The cost functions for the networks we consider in this 
paper make use of more general norms than simply the 
Euclidean norm. Hence, we introduce a generalisation 
of Euclidean spaces, namely finite-dimensional normcd 
spaces or Minkowski spaces. See il2j for an introduction 
to Minkowski geometry. 

A Minkowski space (or finite- dimensional Banach space) 
is K" endowed with a norm ||-||, which is a function ||-|| : 
M" ^ M that satisfies 

• llxll > for aU X e M", = only if x = 0, 

• ||aa;|| = |a|||x|| for all a G M and x G K", and 

• + < M + \\y\\. 

We now discuss some aspects of the SMT problem, since 
this is a special case of the GNP, where all flow are 
zero. Our terminology for the SMT problem is based 
on that used in [5]. Let T be a network interconnecting 
a set N = {pi, . . . ,Pn} of points, called terminals, in a 
Minkowski space. A vertex in T which is not a terminal is 
called a Steiner point. Let G{T) denote the topology of T, 
i.e. G(T) represents the graph structure of T but not the 
embedding of the Steiner points. Then G{T) for a short- 
est network T is necessarily a tree, since if a cycle exists, 
the length of T can be reduced by deleting an edge in the 
cycle. A network with a tree topology is called a tree, its 
links are called edges, and its nodes are called vertices. 
An edge connecting two vertices a, 6 in T is denoted by 
ab, and its length by ||a — 

The splitting of a vertex is the operation of disconnecting 
two edges av,bv from a vertex v and connecting a,b,v 



to a newly created Steiner point. Furthermore, though 
the positions of terminals are fixed, Steiner points can be 
subjected to arbitrarily small movements provided the 
resulting network is still connected. Such movements are 
called perturbations, and are useful for examining whether 
the length of a network is minimal. 

A Steiner tree (ST) is a tree whose length cannot be 
shortened by a small perturbation of its Steiner points, 
even when splitting is allowed. By convexity, an ST is 
a minimum-length tree for its given topology. A Steiner 
minimum tree (SMT) is a shortest tree among all STs. 
For many Minkowski spaces bounds are known for the 
maximum possible degree of a Steiner point in an ST, 
giving useful restrictions on the possible topology of an 
SMT. For example, in Euclidean space of any dimension 
every Steiner point in an ST has degree three. Given a set 
N of terminals, the Steiner problem (or Steiner Minimum 
Tree problem) asks for an SMT spanning N. 

Gilbert [4] proposed the following generalisation of the 
Steiner problem in Euclidean space, which we now extend 
to Minkowski space. Let T be a network interconnecting 
a set N = {pi, . . . ,pn} of n terminals in a Minkowski 
space. For each pair Pi,Pj, i ^ j of terminals, a non- 
negative flow tij — tji is given. The cost of an edge e in 
T is w{te)le, where le is the length of e, is the total 
flow being routed through e, and weight w{-) is a unit 
cost function satisfying 

w{t) > and w{t) > if t > (1) 

W{ti+t2) > w{ti) Vt2>0 (2) 

w{ti+t2) < w{ti)+w{t2) yti,t2>0 (3) 
w{-) is concave (4) 

A network satisfying Conditions ([1]) - ([3]) is called a 
Gilbert network. For a given edge e in T, w{te) is called 
the weight of e, and is also denoted simply by Wg. The 
total cost of a Gilbert network T is the sum of all edge 
costs, i.e. 

G(T) - ^«;(te)Ze 

eeE 

where E is the set of all edges in T. A Gilbert network 
T is a minimum Gilbert network (MGN), if T has the 
minimum cost of all Gilbert networks spanning the same 
point set N , with the same flow demands tij and the same 
cost function w{-). By the arguments of 0|, an MGN 
always exists in a Minkowski space when Conditions ^ 
- (U) are assumed for the weight function. 

Conditions H]) - © above ensure that the weight func- 
tion is non-negative, non-decreasing and triangular, re- 
spectively. These are natural conditions for most appli- 
cations. Unfortunately, the first three conditions alone 
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do not guarantee that a minimum Gilbert network is a 
tree. An example of such a minimum network where the 
flow necessarily splits is given in |16j. (Note that in [3], 
Condition ([3]) , which we call the triangular condition, was 
incorrectly interpreted as concavity of the cost function.) 

The Gilbert network problem (GNP) is to find an MGN 
for a given terminal set N, flows tij and cost func- 
tion w{-). Since its introduction in [3], various aspects 
of the GNP have been studied, although the emphasis 
has been on discovering geometric properties of MGNs 
(see [3], [TT], [T3], [III)- As in the Steiner problem, addi- 
tional vertices can be added to create a Gilbert network 
whose cost is less than would otherwise be possible, and 
these additional points are again called Steiner points. A 
Steiner point s in T is called locally minimal if a per- 
turbation of s does not reduce the cost of T. A Gilbert 
network is called locally minimal if no perturbation of the 
Steiner points reduces the cost of T. 

The special case of the Gilbert model that is of interest 
in this work is when N — {pi , . . . , p„ g} is a set of termi- 
nals in a Minkowski space, where pi, . . . ,p„ are sources 
with respective positive flows ti, . . . , t„, and q is the sink. 
All flows are between the sources and the sink; there are 
no flows between sources. It has been shown in [11] that 
concavity of the weight function implies that an MGN 
of this sort is a tree. Hence we refer to an MGN with 
this flow structure as a minimum Gilbert arborescence 
(MGA), and, as mentioned in the introduction, we re- 
fer to the problem of constructing such an MGA as the 
Gilbert arborescence problem (GAP). 

If vi and V2 are two adjacent vertices in a Gilbert ar- 
borescence, and flow is from vi to V2 then we denote the 
edge connecting the two vertices by viV2. 



deflned as follows: 
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We say that a Minkowski space (R", ||-||) is smooth if the 
norm is differentiable at any x ^ a, i.e., if 

\\x + th\\-\\x\\ 



lim 



t 



exists for all x, /i G R" with x ^ o. It follows easily 
that fx is a linear operator fx : M" — > M and so can be 
represented by a vector x* G R", called the dual vector of 
X, such that {x*,y) = f^iy) for ally £ R", and ||a;*||* = 1. 
In fact X* is just the gradient of the norm at x, i.e., 
X* = V||a;||. 

More generally, even if the norm is not differentiable at 
X, a vector x* G R" is a dual vector of x if x* satis- 
fies {x*,x) — \\x\\ and ||a;*||* — 1. By the Hahn-Banach 
separation theorem, each non-zero vector in a Minkowski 
space has at least one dual vector. A Minkowski space 
is then smooth if and only if each non-zero vector has a 
unique dual vector. 

A norm is strictly convex if ||a;|| = ||?;|| — 1 and x ^ y 
imply that ||^(a; + y)\\ < 1, or equivalently, that the unit 
sphere 

^(||-||) = {2;eR":||x|| = l} 

does not contain any straight line segment. A norm ||-|| 
is smooth [strictly convex] if and only if the dual norm 
II -11* is strictly convex [smooth, respectively]. 

Theorem 1. Suppose a smooth Minkowski space 
(K", [[•[[) is given together with a concave weight function 
w, sources pi, . . . ,p„ G R", and a single sink q £ R", all 
different from the origin a. Let the flow associated with 
Pi be ti. ( See Figure [7J ) For each pi let p* denote its 



3 Characterisation of Steiner Points 

In this section, we generalise a theorem of Lawlor and 
Morgan [B] to give a local characterisation of Steiner 
points in an MGA. The characterisation in [6] holds for 
SMTs, which correspond to the case of MGAs with a con- 
stant weight function. Their theorem is formulated for ar- 
bitrary Minkowski spaces with differentiable norm. Our 
proof is based on the proof of Lawlor and Morgan's the- 
orem given in [8] . A generalisation to non-smooth norms 
is contained in [TO] for SMTs and in [TF for MGAs. Such 
a generalisation is much more complicated and involves 
the use of the subdifferential calculus. 

We first introduce some necessary definitions relating to 
Minkowski geometry, in particular with relation to dual 
spaces. For more details, see [12) . 

We denote the inner product of two vectors x,y € R" 
by {x,y). For any given norm ||-||, the dual norm |)-]|* is 
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Figure 1: A Gilbert network with star topology. 

dual vector, and let q* denote the dual vector of q. Then 
the Gilbert arborescence with edges opi, i = 1, . . . , n and 
oq, where all flows are routed via the Steiner point o, is 
a minimal Gilbert arborescence if and only if 



(5) 



i=l 



and 



[^ w{t,)p* II* < wiY^ t,) for all I Q {!,..., n}. (6) 
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Note: We think of Condition [5] as a flow-balancing condi- 
tion at the Steiner point, and Condition [6] as a condition 
that ensures that the Steiner point does not spht. 



Proof. We are given that the star is not more ex- 
pensive than any other Gilbert network with the same 
sources, sink, flows and weight function. 

In particular, o is the so-called weighted Fermat-Torricelli 
point of pi, . . . ,p„, g with weights ii, . . . , t„, X^Li re- 
spectively, which implies the balancing condition ([5]) . We 
include a self-contained proof for completeness. If the 
Steiner point o is moved to —te, where t gM. and e € R" 
is a unit vector (in the norm), the resulting arborescence 
is not better, by the assumption of minimality. Therefore, 
the function 

n 

n 

+ w{Y^U)(\\q + te\\-\\q\\)>Q 
1=1 

attains its minimum aXt~Q. For t in a sufficiently small 
neighbourhood oi Q, pi + te ^ o and q + te ^ hence tp^ 
is differentiable. Therefore, 



\\P^ 



te\ 



= v'l(0) = liin ^ 

\\q + te\\ 



W 



i=l 



= E^(i^)(p:,e)+u;(EtO {q*,e) 

1=1 i=l 

(n n \ 

1=1 i=l I 

Since this holds for all unit vectors e, ([5]) follows. 

To show ([SJ for each / C {!,..., n}, we may assume 
without loss of generality that 1^0 and / 7^ {1, . . . , n}. 
Consider the Gilbert network obtained by splitting the 
Steiner point into two points o and +te € R, e a unit 
vector) as follows. Each pi, i ^ /, is still adjacent to o 
with flow ti, and q is joined to with flow X]r=i 
now each p^, i G /, is adjacent to te with flow t^, and te 
is adjacent to o with flow J2iei^i^ shown in Figure [S] 

Since the new network cannot be better than the original 
star, we obtain that for any unit vector e, the function 

n 



iei 



i=l 




Pi,i ii 
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Figure 2: The Gilbert network obtained by splitting the 
Steiner point o. 



ferentiable at 0, we can still calculate as follows: 



< lim 

t^o+ t 



lim > w(ti) 



lb* " - \pi 



w 



j=i 



= {Y.w{t,)p*,-e)+w{J2t^)■ 



1=1 



Therefore, (I]ie/ '"(^Op*- e) < ^(X^ILi ^0 ^^li* 
vectors e, and (El follows from the definition of the dual 



norm. 



(-4=) Now assume that pl,Pn,Q are dual unit vectors that 
satisfy ^ and ([5]). Consider an arbitrary Gilbert ar- 
borescence T for the given data. For each i, let Pi be 



the path in T from pi to (7, i.e.. Pi = x^'^X2'^ 



(i) 



where x^ ^ = Pi,x')^^ = q, and ^^''^^j+i are distinct edges 
of T for j = 1, . . . , fci — 1. For each edge e of T, let 
S'e = {i : e is on path Pi}. Then the flow on e is X^iGS ^« 
and the total cost of T is 



E ^"(E^')ii^-2^ii 



e—xy is iGSe 
an edge of T 



The cost of the star is 



j2wit,)\\p,\\+w{Y,tiM\ 

n n 

= ^(^0 {plP'd + wiY, {q*,q) 

1=1 i=l 
n 

= 5]]u;(t,)(p:,P<-'Z) by® 



attained its minimum at t = 0. Although ip^, is not dif- 



1=1 

n fci — 1 

E^(^^)E(?'*' 



1=1 



4 



e—xy is \i^Se I 
an edge of T 

e—xy is i^Se 
an edge of T 

< E «^(E*^)ii^-yii by®- □ 

e—xy is i^S^ 
an edge of T 

Note that the necessity of the conditions ([5]) and ^ holds 
even if the weight function is not concave. It is only in the 
proof of the sufficiency that we need all minimal Gilbert 
networks with a single source to be arborescences. 

4 Degree of Steiner Points in a 
Minkowski plane with linear weight 
function 

We now apply the characterisation of the previous section 
in the two-dimensional case, assuming further that the 
weight function is linear: w{t) ~ d + ht, d > 0, h > 0. 

Theorem 2. In a smooth Minkowski plane and assuming 
a linear weight function w{t) = d + ht, d > 0, h > 0, a 
Steiner point in an MGA necessarily has degree 3. 




Figure 3: A polygon with edges corresponding to v*. 




Pn 



Figure 4: Illustration of the proof of Theorem 2 for a unit 
ball with straight line segments on the boundary. 



Proof. By Theorem [TJ an MGA with a Steiner point of 
degree n + 1 exists in with a smooth norm ||-|| if and 
only if there exist dual unit vectors pj, . . . ,p* , g* G 
such that 

n n 

^(d + ht,)p* + {d + hY U)q* - o 

i=l i=l 

and 

\\J2id+hU)p*\\* <d+hJ2t^ for all /C {!,...,«}. 

is/ iei 

Label the p* so that they are in order around the dual unit 
circle. Let v* = {d + ht^)p* and w* = {d + h^,"^^ t,)q*. 
Then the conditions become 

vl + ■ ■ ■ + + w* — o, 

and 

^ '^ + ^Y*^ for all / C {!,..., n}. (7) 

Thus we may think of the vectors u^f, . . . , u* , w* as the 
edges of a convex polygon with vertices a* = J^l^iV*, 
J = 0, . . . , n in this order (see Figure [3]). 

Assume for the purpose of finding a contradiction that 
n > 3. Then the polygon has at least 4 sides. Note that 
the diagonals QqU* and a*_]^a* intersect. Applying the 



triangle inequality to the two triangles formed by these 
diagonals and the two edges v* and w* (as illustrated in 
Figure [3]), we obtain 

||a;r + ||a:-a*_jr>||t-*|r + I|u;*r 

n 

^ d + htj + d + h^ti 

i=l 

3 n 

= d + hy^Ji + d + hy^^tj 

1=1 i=j 
3 n 

>iiE<ir + iiE<ir ^y© 

1=1 i=j 

= \\a*r + \\al~a*_,r. 

Therefore, equality holds throughout, and we obtain 
equality in the triangle inequality. Since the dual norm 
is strictly convex, it follows that v* and w* are paral- 
lel. This holds for all j = 2, . . . , n — 1. It follows that 
p| = • ■ • = p* = —q*. Geometrically this means that the 
unit vectors jj^Pi and ~ 111119 all have the same support- 
ing line on the unit ball. We can think of this condition 
on the vectors pi and g as a generalisation of collinearity 
to Minkowski space. 

Choose a point S2 on the edge opi such that the line 
through S2 parallel to op„ intersects the edge opi in si, 
say, with si ^ a. See Figured! Because of the straight 
line segments on the boundary of the unit ball, ||a;-f t/|| = 
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1 1 a; 1 1 + [|y|| for any x,y such that the unit vectors 
and ji^y he on this segment. In particular, 

||S2 - Sill + oil = 11^2 -0||. (8) 

Now replace p^o by the edges P2S2 and S2S1 , replace pio 
by pisi and sio, and add the flow t2 to sio. The change 
in cost in the new Gilbert arborescence is 

{W{tl)\\pi -Sl\\+ W{tl + <2)||S1 - o|| 

+ W{t2)\\p2 - S2II + W{t2)\\s2 - Sill) 
- {IV{U)\\PI - o\\ - W{t2)\\p2 - 0\\) 

= -w;(ti)||si|| - W(t2)||s2|| +t2)||si|| 
+ z«(t2)(||s2||-||si||) by® 

= {w{ti + t2) - W{ti) - W{t2))\\si\\ 

= {d + h{ti+t2) - (d + hti) - (d + M2))||si|| 
= -d||si|| <0. 

We have shown that a Steiner point of degree at least 4 
leads to a decrease in cost. In an MGA a Steiner point 
must then necessarily be of degree 3. □ 

5 Conclusion 

In this paper we have studied the problem of designing 
a minimum cost flow network interconnecting n sources 
and a single sink, each with known locations and flows, 
in general finite-dimensional normed spaces. The net- 
work may contain other unprescribed nodes, known as 
Steiner points. For concave increasing cost functions, a 
minimum cost network of this sort has a tree topology, 
and hence can be called a Minimum Gilbert Arbores- 
cence (MGA). We have characterised the local topologi- 
cal structure of Steiner points in MGAs for linear weight 
functions, specifically showing that Steiner points neces- 
sarily have degree 3. 
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